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Analysis of the two-point velocity correlations 
in turbulent boundary layer flows 

By M. Oberlack 
1. Motivation and objectives 

Two-point Rapid Distortion Theory (RDT) has become an important tool in the 
theory of homogeneous turbulence. Modelers try to implement appropriate results 
from RDT in their statistical turbulence models, for example in the structure based 
model developed by Kassinos and Reynolds (1994). 

On the other hand, in non-homogeneous equilibrium flows the logarithmic law 
is one of the cornerstones in statistical turbulence theory. Experimentalists have 
found the log-law in a broad variety of different turbulent wall shear flows, and 
statistical models have been made to be consistent with the log-law. 

The logarithmic law was first derived by von Karman (1930a, 1930b) using di- 
mensional arguments. Later Millikan (1939) derived the law-of-the-wall more for- 
mally using the so called “velocity defect law”, also first introduced by von Karman 
(1930b). Even though the derivation was much more comprehensive from a phys- 
ical point of view, the velocity defect law is essentially an empirical observation. 
A first derivation of the law-of-the-wall using asymptotic methods in the Navier- 
Stokes equations was given by Mellor (1972). Mellor needed the viscous sub-layer 
to obtain the log-region, and his scaling of the inertial range in the log-region is in 
error because it does not give the one-point limit of production equals dissipation. 

The general objective of the present work is to explore the use of RDT in analysis 
of the two-point statistics of the log-layer. RDT is applicable only to unsteady 
flows where the non-linear turbulence-turbulence interaction can be neglected in 
comparison to linear turbulence-mean interactions. Here we propose to use RDT 
to examine the structure of the large energy-containing scales and their interaction 
with the mean flow in the log-region. 

The contents of the work are twofold: First, two- point analysis methods will 
be used to derive the law-of-the-wall for the special case of zero mean pressure 
gradient. The basic assumptions needed axe one-dimensionality in the mean flow 
and homogeneity of the fluctuations. It will be shown that a formal solution of the 
two-point correlation equation can be obtained as a power series in the von Karman 
constant, known to be on the order of 0.4. 

In the second part, a detailed analysis of the two- point correlation function in the 
log-layer will be given. The fundamental set of equations and a functional relation 
for the two-point correlation function will be derived. An asymptotic expansion 
procedure will be used in the log-layer to match Kolmogorov’s universal range and 
the one-point correlations to the inviscid outer region valid for large correlation 
distances. 
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2. Governing equations of the two-point velocity correlation function 

Using the standard Reynolds decomposition U, = «, + u, and P = p + p, the 
Reynolds averaged Navier Stokes (RANS) equations read 

dui _ dui 1 dp (Put dujui n , 

dt +Uk dxk pdii V dx\ dxt 

and the fluctuation equation, later on referred to as Adequation, is 


du< dui dui 

Mx)= — +ut — + u t — 


dujuT dujUk 1 dp 

dxk dxk p dii 


d 2 u, _ 
V dxl 


The corresponding continuity equations are 


du k 


= 0 and 


du k 


= 0 . 


( 2 ) 


( 3 ) 


dxk ~ ~~ dxk 

The five two-point correlation tensor functions that appear in the two-point cor 
relation equation (5), further below, are defined as 


R ij (x 1 r\t) = u i (x,t)u j (xW,t) , 
puj(x,r;t) = pOM)“j(* (1) >*) » 
ujp(x,r;t) - u j (x,t)p(x( 1 \t) , 

R(ik )} (x,r-,t) = Ui(x,t)u k (x,t)uj(xW ,t) , 

R i{jk) (x,r;t) = Ui(x,t)uj(xW,t)uk(xW,t) . (4) 

All tensors in (4) are functions of the physical and the correlation space coordinates 
x and r = x^ — x respectively. The double two-point correlation Rij, later on 
simply referred to as two- point correlation, converges to the Reynolds stress tensor 
ujuj in the limit of zero separation r. 

The well known two-point correlation equation (Rotta (1972)) can be written as 


DRij „ diii(x,t) 
~dT ~ Kkj dx k 


— Rik 


dujjXjJ) 

dx k 


i x+r 


[u fc (* + r,<)-u*(*,t)] 


dRy 

dr k 



dpuj dujp 
dri + drj 


+ v 


d 2 Rij 2 d 2 Rj } cPRjj ' 

dxkdxk dxkdr k dr k dr k _ 


21 , 2 

dx k dr k 


[•R(it)> - R i(jk)} 


( 5 ) 


For both two-point velocity- pressure correlations, ujp and puj a Poisson equation 
can be derived. The divergence d/dii - d/dr, of equation (5) leads to a Poisson 
equation for puj, 
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Figure 1. 


Sketch of the coordinate system and the mean velocity field 


1 

P 


gVj 2 &puj 

dxkdxk drkdxk 

dxkdxi dxkdri 


d 2 pu } ' 
dr k dr k 
d 2 Rjki)j 
dr k dr i 


2 


_ dR l} ' 

dx/ [ c?x* 5rjt 


( 6 ) 


and the divergence d/drj leads to the corresponding Poisson equation for Uip, 


i d 2 u, P _ 2 a 2 i?, (fc/) 

pdr k dr k dxi x+r dr k dr k dr t 


(7) 


where the vertical line means that the derivative is taken with respect to x but will 
be evaluated at x+r. All of the dependent variables in (5)— (7) have to satisfy the 
continuity conditions 


dR,j dRij 

dx{ dr, 



and 


0 , 

*1 

ii 

o 

(8) 

duip 

dv~p 

(9) 

dxi 

*7 = 0 - 


For the analysis of the self-similar, two- point correlation equation further be- 
low, two identities are important. They can easily be derived from a geometrical 
consideration by interchanging the two points * and x (1) = * + r 


Rij(x,r) — Rj,(x + r,-r) , Uip(x,r) = pui(x + r,— r) . (10) 

The latter identities are the key elements for the derivation of some boundary con- 
ditions and for the deeper understanding of the self-similar two-point correlations. 
There exists a similar identity for the triple correlation which will not be used here. 
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3. The log-law - a self-similar form of the two-point correlation equation 

A sketch of the coordinate system and the mean velocity field adopted in the 
proceeding paper is given in Fig. 1. Within this subsection it will be shown that 
the logarithmic part of the law-of- the- wall mean velocity profile can be derived 
from the two-point correlation equation and hence from the Navier- Stokes equation 
if there exists a regime where the following assumptions hold: 

• the mean velocity is parallel to the wall; 

• the statistics in that domain are independent of viscosity and time; 

• the Reynolds number is high; 

• no mean pressure acts on the flow field. 

The last assumption can be eliminated, but in this approach it will be focused on 
the zero pressure gradient case. Beside the above assumptions no other conditions 
are needed in order to determine the log-law mean velocity profile and the self- 
similarity of the correlation functions. 

Inferring the above assumption in the Reynolds averaged Navier-Stokes equations 
(1), it is easy to confirm that the gradient of the Reynolds stress tensor on the right 
hand side is the only remaining term. Integrated one time we obtain that Hi u] is 
independent of x. However, this is not necessarily true for the two-point correlation 
tensor R xy It could depend on x if the dependence vanishes in the zero separation 
limit. This can only be achieved by having the following dependence on a new 
variable 


r = rg(x) (11) 

where g(x) has to be determined later and no other hidden dependence on x can be 
in the correlation functions. Of course, the latter definition of f can be generalized 
to different unknown scaling functions for every component of r, but from equation 
(5) it can be verified that only a single scaling function exists. With the above given 
assumptions, defining x = x and using the transformation rules 

d d 1 dg „ d 

dx t dxi + g dii Tk dr* 

the i?,j-e quation (5) reduces to 


A- A 

dr k 9 df k 


( 12 ) 


0 = 

ax-2 


dx 9 


x 2 +r 2 


, , N - / M d Ri i 

- [«i (x 2 + r 2 ) ~ u i (i 2 )J 9 


1 dg _ dpuj 
P [fif dii k dfk 


dpu] duip 

9 -^7— + 9 
9r, OTj 


1 dg „ dR(ik)j , ^ d D 1 

gdTk n ~dfT + 9 dfk [ R(tk)j R,Uk " 


( 13 ) 
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As mentioned above, there is no hidden x dependence in the correlation function 
and therefore all x derivatives coming from (12a) have been omitted. Obviously, 
equation (13) can only have a non-trivial solution, and thus be independent of X, if 
all the coefficients have the same functional dependence on x. Hence, the following 
set of differential equations determine the u \ and g dependence on * 


and 



for t = 1 , 2,3 , 


<fal(*2) 

dx 2 


~ 9 


(14) 


1 

9 


du i(x 2 ) 


dx 2 


x+r 


= /i (*)*/a(«) 


u i( x 2 + r 2 ) - ui(x 2 ) = f 3 (f) ^ fi(x) 


(15) 


are additional consistency conditions for Uj and g. The last equation in (14) de- 
termines g to depend only on x 2 . Hence, the equations (14) have two independent 
sets of solutions given by 


9(x 2 ) = 


(i)/ 1 

c 2 (l 2 - c\ 


< 1 )> 


and 


JD 

M x i) = -pyln(* 2 -<4 1) ) + ci 1> 


(16) 


< 7 (^ 2 ) = 4 2) , tll(* 2 ) = 4 2 >c S 2 ) *2+4 2) (17) 

where the c^ p) ’s are integration constants or proportionality factors. Obviously, only 
the first set of equations correspond to a boundary layer flow because the solutions 
(17) define homogeneous shear turbulence which contradicts the assumption to be 
independent of time. Both equations (15) require c( J) = 0 and can be absorbed 
in the correlation functions. In common notation we finally obtain 


and 


Hi — u r 6 tl 


— ln(i 2 ) + C 


(18) 


where u T is defined as 



(19) 


U r = 



(20) 


Inserting (18) and (19) into equation (13) and multiplying by the von Karman 
constant k the final form of the /£,y -equation results: 
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where 


0 — ^*1 ^ i 2 


1 


K 

+ - 
p 


+ K 


1+^2 

dpu] 


dR*i 

- ln(l + r 2 ) -^r 1 
Or i 


s„rM :+ 

arfc ar 

- ^(,2)> . g 

1 dft dfk 


r_^ni 

i af> J 

^ (*?« M - ***>)] 


( 21 ) 


/?*, = ^ 






ft 




R i(jk ) ~ 


R 


'»(>*) 


. P u ) 


j * 

and Ujp = — v 


u*P 

.3 


( 22 ) 


The procedure described above can be extended to the three-point triple-corre- 
lation equation said any higher order correlation equation if an additional spatial 
point is introduced for each additional tensor order. As a result it is easy to verify 
that the whole set of equations define an infinite set of linear tensor equations but 
which are far too complex to be solved in general. Nevertheless, it is worthwhile to 
analyze some features of the solution. 

In principle this infinite set of equations could be solved by the following pro- 
cedure. Beginning with the two-point correlation equation, the triple correlation 
can be considered as an inhomogeneous part of the Rij equation. Once the ho- 
mogeneous solution is obtained, the inhomogeneous solution can be computed by 
standard methods. In the next step, the triple- correlation equation has to be tack- 
led and its solution will be substituted in the solution for and so forth for 
higher correlations. In each equation the von Karman constant k only appears as 
a factor of the highest order tensor and hence the final solution for Rij is a power 
series in k 


r-, = E 4"’““ • < 23 > 

m=0 

represents the solution of the two-point correlation equation after neglecting 
the triple-correlations and all higher order terms. 

The structure of the formal solution in equation (23) admits the hope that a 
truncated series may provide some insight in the log-law statistics. Hence, in the 
following the triple-correlations will be neglected. Using the similarity variable in 
the poisson and the continuity equations, the puj equation (6) becomes 


. . a 2 puj- cp-pui' . Pf*]- ap»r „«i=r 

k ' avT + ~dhh + 2 ‘ avT + ‘"wT + 


2 p dR* 2] 

k df\ ’ 


(24) 


the Uip equation (7) becomes 
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d 2 uip* 2 p i dm 


>2 


df k f k k 1 + r 2 dr i 
and the continuity equations (8) and (9) yield 

dR*, 


. dm 2] 

r k^r L + 
dr k 


^=0 , 


dr, 


dR*, 


dpu 7 n , - du5p 

o- = 0 and r k -—— 
dr j dr k 


+ 


df i 

dujp* 

dr. 


= 0 


= 0 


(25) 


(26) 

(27) 


The identities (10) can also be transformed in a similar manner. Introducing the 
transformation (19) into the equation (10a), we obtain the relation Rij(x 2 ,x 2 r) = 
■Rjii(z 2 (l + r 2 ), -x 2 r). Because it was previously assumed that rill two-point corre- 
lation functions are solely functions of r, only the ratio of the first and the second 
parameter can appear in R tJ . This argumentation can be extended to the pressure 
velocity correlation. Thus, we finally obtain 


(r+y < 28 > 

and 

u^*(f) = pu 7* ■ (29) 

The latter identity also holds if xTip m and pu , * are interchanged. 

These two relations give valuable insight into the structure of the solution. Rela- 
tion (28) connects different f domains to each other and provides boundary condi- 
tions in the f 2 direction. 

One interesting feature of (28) is that it can be considered as a functional equation 
for each trace element. It is easy to verify that one solution, but probably not the 
most general solution to equation (28), is given by the following form 


*[Vr) ( f ) = F -r 


(m^>)’.M 


r 2 r 2 


(30) 


where is one of the three trace elements of R* } . 

In addition if the solution for any off-diagonal R* } element (i ^ j) is known, (28) 
provides the solution for the R* t . A similar feature for u~p t and pu 7 is given by 
relation (29). 

If boundary conditions have to be satisfied in infinity, all correlation functions 
decay to zero. Therefore, any solution of equations (21) and (24)-(27) have to obey 
the boundary conditions 


R*iji f k — ► ±oo) = pui*(r k -v ±oo) = uip*(r k — ► ±oo) = 0 for k = 1,3 (31) 
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x 2 

x 2 = — r 2 



FIGURE 2. Sketch of the boundary condition in the x 2 -r 2 plane, 
and 


R*j(f 2 -*■ oo) = puj*(?2 -*■ oo) = u7p*(f 2 -» oo) = 0 . (32) 

To better understand the boundary conditions in the wall-normal direction, a 
sketch of the x 2 -r 2 plane is given in Fig. 2. Picking any value for x 2 , the negative 
part of r 2 can not be smaller than x 2 and hence one bound is on the line x 2 = — r 2 . 
The bound for the physical coordinate is at x 2 = 0. Using the definition of the 
scaled non-dimensional coordinate (19), it is clear from Fig. 2 that r 2 represents the 
inverse of the slope given by any straight line through the origin ranging between 
the two latter bounds. Hence, the domain for f 2 is restricted to — 1 < r 2 < oo. 
Using (28) and (29) together with (32) one obtains 

R*j( ?2 = -1) = o (33) 


and 


pu*(f 2 = -1) = u,p*(f 2 = -1) = 0 . (34) 

Obviously, the boundary conditions are all homogeneous and one may expect the 
solution to be zero. In section (5) it will be discussed why the equations might have 
a noil-trivial solution, but a rigorous proof is still outstanding. In the next section 
an integral relation coming from the one-point equations will be derived, which 
closes the missing information regarding the scaling of the two-point correlations. 

4. Kolmogorov’s universal range and one-point correlations 

The self- similarity of the correlation functions introduced in section 3 is only valid 
in the limit of large Reynolds number, based on the wall distance and the friction 
velocity 


Rt T = 


u r x 2 


V 


(35) 
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This is also the definition of y + . From experiments it is known that the log region 
starts at about y + = 40 and extends to y + = 0.2U6/v. 

The analysis in the previous chapters is inviscid, and hence is not a regular 
expansion in Re r . It is not applicable for small correlation distances, as will be 
explained in some detail now. An inner viscous layer in correlation space has to 
be introduced in order to meet the requirement that viscosity is important for the 
dissipation tensor £,j in the one-point limit. 

Comparing the two-point correlation equation (5) in its most general form to the 
inviscid version in the log-layer (21), no viscous term has been retained. In contrast 
to that, the Reynolds stress transport equation in the log-layer 

- [uiU 2 6ji + uTujiii] — + $ij - £ij = 0 (36) 

KX'i 

contains the viscosity v in the dissipation tensor, defined by 


€ij 


= 2v 


dui duj 
dx k dx k 


= 2v lim 
r=0 


d 2 R,j 

dx k dr k 


and the pressure-strain tensor is defined by 


PR,; ' 

dr k dr k 


(37) 


p 

duj du t 

= lim 

dpitj du x p dujp 

p 

dx x + dxj 

r=0 

dr x dxj drj 


(38) 


The contraction of equation (36) together with U 1 U 2 = —u\ determines the scalar 
dissipation 


£ = 


£hk _ 

2 KX 2 


(39) 


As mentioned above we find from equation (36) that the asymptotic arguments 
we have used so far are not valid for correlation distances on the order of the 
Kolmogorov length scale l v . The Kolmogorov length and velocity scale are given by 


U = j = x 2 Re T * and u v = = u r Re r 7 k J . (40) 

The only scaling of the independent variables with which the correct balance can 
be achieved in the two-point correlation equation is given by 


r Sir 

£ = 7- = Re* k i — 

*n Xo 


(41) 
X2 

In line of Kolmogorov’s arguments, the seeding of the dependent variables must be 


Rij = ulw]' - < [4 0) (O + o (ife;*)] , 
R(tk)j = u r, D(i k )j(£) + O (i?e r , 
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Ri{jk) — u n [ 
pu] = pu\ 
Ujp = pu 3 v 


<»,(«) + o(ju;»)] , 
[m} 0) («) + o (&;*)] , 
[jv,< 0 '«)+o (&;*)] . 


( 42 ) 


Putting (41) and (42) into (5), (6) and (7) the leading order terms in each equation 
are given by 


dM { ° ] 

-S,iu 2 uj* - 6jiUiU 2 * + — 7j7 — 

* 


ON) 


(o) 


d 2 B\f 
-2 1 


dD 


+ 


(o) 

(»*)> 


dD K 


(o) 


i(jk) 


d(,j di k di k dik 


dtk 


(43) 


and 


d£kd(,k d(, k dii 


d 2 Nl 0) = d 2 D ( °l t) 
d£kd£k d£k d£i 


(45) 


In order to obtain a uniform solution there must an overlapping region that 
matches the inner and the outer solution together. From (42a) we see that the limit 
£ — ► oo in the inner layer of the two-point correlation converges to the Reynolds 
stress tensor and the same must be valid for a solution of the equations (21) and 
(24)-(29) in the outer layer for the limit f —* 0. Using the same limits for both 
regions in the triple- and the pressure- velocity correlations, they both drop to zero 
as they should do. As a result, the matching between the inertial subrange and 
obviously specifies the outer solution at r = 0 to be but the actual 

numerical value of Reynolds stress tensor is still unknown. 

Note, that the equation corresponding to (43) in Mellor’s paper (1972) (his equa- 
tion (59)) has a serious error. It does not have the production terms which, of 
course, axe responsible for the energy transfer rate. 

As mentioned above the inner layer does not determine the absolute value of 
the Reynolds stress tensor because the triple correlations can not be neglected in 
(43)-(45). Thus an additional assumption is needed to determine the values of 
In Kolmogorov’s original hypotheses it was suggested that in the limit of large 
Reynolds number the dissipation will be isotropic. Saddoughi’s (1994) very high 
Reynolds number experiment of a turbulent boundary layer in a wind tunnel sup- 
ports this idea of isotropy. Hence, we take 


£ij = ^ Sije . (46) 

Using this, the three trace elements of 4> u can be obtained from the Reynolds Stress 
tensor equation which in non-dimensional form can be written as 
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- [nTuT <5,1 + + $*, - = 0 with **> = *ij^ (47) 

O U T 

or in component notation 

*11 = , *2*2 = f > * 3*3 = l • (48) 

Note, that the latter result for the pressure-strain correlation holds no matter what 
is assumed for the triple-correlations. As a result, all high Reynolds number second- 
moment closure models should be consistent with this result. In most second mo- 
ment models this could only be ensured by adding wall reflection terms to the 
pressure-strain model. 

Because the system (21) and (24)-(29) has homogeneous boundary conditions 
on all boundaries, there is nothing that specifies the amplitude of R*j or the value 
of UiUj* as mentioned above. In fact, this would also be true if higher correlation 
functions would have been taken into account. The definition (38) together with 
the result (48) can be used to calculate the values for the Reynolds stress tensor. 

The term on the right-hand side of (38) can be rewritten as an integral of the 
two-point correlation and some boundary integrals. This was necessary because the 
limit r — > 0 has to be evaluated within the dissipation range where not enough is 
known about the two-point velocity-pressure correlation. It can be found that the 
dissipation range, which is of the order of Z^, makes a higher order contribution to 
the above mentioned integral in the limit of large Reynolds number and thus can 
be neglected. After neglecting the triple-correlations we find 


<r 



i 

1 + r 2 



+ n 


d 2 R'u \ 

dr\d?\ ) 


dfjdf\\ |r| 


J) (49) 


where (i +-» j) abbreviates the addition of the previous term with indices inter- 
changed. No boundary integral has to be kept due to the homogeneous boundary 
conditions for all variables. Once a solution to the equations (21) and (24)-(27) are 
computed the scaling of the two-point correlations can be calculated by equating 
(48) and (49). Using this, the value for u t u]* can be taken from at r = 0 as 
has been proven by the matching between the Kolmogorov universal range and the 
outer inviscid solution. 


5. Future plans 

There are basically two outstanding problems within the whole approach of RDT 
in the log-layer. The first one is the fact that it has to be proven that the system 
(21), (24)-(29) has a non-zero solution even though all boundary conditions are 
homogeneous. A strong hint towards this character of the equation is gained by 
the analysis of the discretised set of equations which, of course, is linear. To see 
why the equations may have a non-zero solution, a result from linear algebra may 
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be recalled. If in a linear system of the form Ax = 0 the matrix A has the rank 
£ and Q < n where n is the number of equations, then the system has nontrivial 
solutions. In this particular case considering the discretised equations (21), (24)- 

(27) , A is a quadratic matrix and its rank can only be smaller than n if there is 
some redundancy in the equations. In fact, this redundancy is due to the identities 

(28) and (29). Even though the structure of the discretised system provides some 
information, the proof of a corresponding feature in the differential equations is still 
outstanding. Once the previous problem is solved, a numerical algorithm has to be 
coded to solve the discretised equations (21) and (24)-(29) because it is very unlikely 
that an analytical solution can be found. In the next step of post-processing the 
numerical results, the ability of the asymptotic limits used in the RDT of the log- 
layer has to be revised and if necessary enhanced by including higher correlations 
in the analysis. Finally, the results of the theory will be compared with DNS data 
from the turbulent channel flow (Kim et al. 1987). 
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